


____ DATASET __________________________________________________________________________________

In this dataset we include 14 files containing the code and data associated to our manuscript
"A Biomathematical Model of Tumor Response to Radioimmunotherapy with aPDL1 and aCTLA4". The
model is coded in MATLAB. The dataset includes 5 .m functions (CODE), 6 .mat workspaces
(DATA) and 1 .ods document containing tables with the experimental data obtained from the
articles by Dewan et al. and Deng et al.

____ MAIN CHANGES FROM V2 _____________________________________________________________________

The immune-death term in the model has been changed to a more consistent one.
A couple of bugs were found and corrected in the implementation.
Complete re-analysis of the experimental data. 


____ CODE _____________________________________________________________________________________

radioinmuno_response_model.m
  Numerical code to solve the radioimmunotherapy response model.
    
annealing_optimization.m
  Simulated annealing to fit model to experimental data (Dewan et al. & Deng et al.).
    
calculate_TCP.m
  Function to calculate TCP for different treatments combining RT & IT.
    
get_equivalentBEDtreatment.m
  Function to calculate treatment schedules with the same BED than a given one.
    
sensitivity_analysis.m
  Function to do the sensitivity analysis.
  
---> Detailed information can be found in the header of those files.


____ DATA _____________________________________________________________________________________

The file name indicates the figure in the paper that was constructed with each workspace data.

Figure_2.mat
  Data for fits to experimental data reported by Dewan et al.: experimental data extracted from
  the article, input parameters for our model, and outputs of our model.
  RT damage model: modified LQ
  Vascular damage: no
  
Figure_3.mat
  Data for fits to experimental data reported by Dewan et al.: experimental data extracted from
  the article, input parameters for our model, and outputs of our model.
  RT damage model: LQL
  Vascular damage: no
  
Figure_4.mat
  Data for fits to experimental data reported by Deng et al.: experimental data extracted from
  the article, input parameters for our model, and outputs of our model.
  RT damage model: modified LQ
  Vascular damage: no

Figure_5.mat
  Data for fits to experimental data reported by Dewan et al.: experimental data extracted from
  the article, input parameters for our model, and outputs of our model.
  RT damage model: LQ
  Vascular damage: yes
  
Figure_6.mat
  Results for the optimization of the treatment with RT and IT in different schedules.

Figure_7.mat
  Results for the investigation of the optimal radiotherapy fractionation.
